Difference in multiplicity distributions in proton-proton and proton-antiproton 

collisions at high energies 
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Secondary charged hadrons multiplicity distributions in proton-proton and proton- 
antiproton collisions differ on principle. There are three types of inelastic processes in 
proton-antiproton scattering. The first type is production of secondary hadrons shower 
at gluon string decay. The second type is shower produced from two quark strings decay, the 
third type is shower produced from three quark strings decay. At the same time there are 
only two types of inelastic processes for proton-proton scattering - gluon string shower and 
two quark strings shower. Theoretical description of multiplicity distributions is obtained 
for proton-proton collisions at energies from 44.5 GeV to 200 GeV and for proton-antiproton 
collisions at energies from 200 GeV to 1800 GeV. The difference between proton-proton 
and proton-antiproton multiplicity distributions is discussed. The predictions of multiplicity 
distribution and mean multiplicity at LHC energy are given. 

PACS numbers: 12.40Nn, 13.85.Hd, 13.85Lg 
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I. INTRODUCTION 

Multiplicity distributions and total cross 
sections are the main characteristics of hadrons 
multiple production at high energies. Analy- 
sis of secondary hadrons multiplicity distribu- 
tions is very important for verification of different 
phenomenological approaches and models. Ex- 
perimental data on multiplicity distribution are 
highly informative and can be measured with suf- 
ficient accuracy which gives significant advantage 
when comparing with theoretical calculations. 



II. LOW CONSTITUENTS NUMBER 
MODEL AND TOTAL CROSS SECTIONS 

We are based on the Low Constituents 
Number Model (LCNM) which was proposed by 
V.A. Abramovsky and O.V. Kancheli in 1980 [if. 
Three steps are concretized in this model: prepa- 
ration of colliding hadrons initial state; interac- 
tion; moving apart of reaction products. 

1) On the first step before the collision there 
is small number of constituents in hadrons. They 
are valence quarks and few gluons which fill in 
the whole spectrum in rapidity space. 



2) On the second step the hadrons interac- 
tion is carried out by gluon exchange between the 
valence quarks and initial gluons and the hadrons 
gain the color charge. 

3) On the third step after interaction the col- 
ored hadrons move apart and when the distance 
between them becomes larger than the confine- 
ment radius, the lines of color electric field gather 
into the string. This string breaks out into sec- 
ondary hadrons. 

In the model total cross sections of proton- 
proton and proton-antiproton collisions are de- 
scribed by the following formulae 



ogp = 63.52s- - 358 T 35.43s- - 56 + 



+<7 + cri Ins + cr 2 (lns) 2 + 
-o"3(lns) 3 + [higher powers of Ins]. 



:1) 
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where first two terms present non vacuum contri- 
butions and were taken from the paper by Cudell 
et. al. 0. There are also constant term and 
terms proportional to powers of logarithm of full 
energy. Constant term corresponds to gluon ex- 
change between initial states containing only va- 
lence quarks. Powers of logarithm are equal to 
number of gluons in initial states of hadrons. 

The fitting of formulae (pQ) to total cross sec- 
tions of pp and pp collisions (Fig.l, data were 
taken from [3]) showed that there are only two 
gluons in initial state. Contribution of the third 
gluon is about 1 mb at energy 14 TeV and thus 
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FIG. 1: Fitting of total cross sections of pp and pp 
collisions in LCNM 



is negligible. The values of parameters are <7o = 
20.08 ± 0.42, o\ = 1.14 ± 0.13, er 2 = 0.16 ± 0.01. 

Having limited number of gluons we can de- 
scribe the types of inelastic processes in proton- 
proton and proton-antiproton interactions. 



III. THREE TYPES OF INELASTIC 
PROCESSES 

From the obtained result of presence of only 
two gluons in initial state it follows existence of 
three types of inelastic processes in pp and pp 
collisions. 

There are three types of inelastic processes 
in proton-antiproton interaction. The first type is 
production of hadrons shower from decay of gluon 
string, it corresponds to constant contribution to 
total cross sections (Fig. 2a). The second type is 
shower produced from decay of two quark strings, 
it corresponds to contributions from one and two 
gluons (Fig. 2b). The third type is shower pro- 
duced from decay of three quark strings, it cor- 
responds to part of contribution from two gluons 
(Fig. 2c). In this case quark strings are produced 
between every quark of proton and antiquark of 
antiproton. 

At the same time there are only two types of 
inelastic processes in proton-proton interaction, 
they are shower from gluon string (Fig. 2d) and 
shower from two quark strings (Fig. 2e,f). There 
are no three quark strings in pp interaction be- 






FIG. 2: Types of inelastic processes for pp (a - c) 
and pp interaction (d - f). Solid lines correspond to 
quarks and antiquarks, wavy lines to gluons, spirals to 
strings. Interaction in final state is marked by dashed 
block. 




a 



FIG. 3: Hadrons production in gluon string, a) Am- 
plitude of transition of two octet states produced after 
gluon exchange to hadrons, dashed lines show direc- 
tion of octet states, b) One of possible diagrams of 
Fig. 2a; dashed block shows one of possible diagram 
of amplitude of transition of two octet states of Fig. 2a 
to hadrons 



cause in this case strings can be formed between 
quark of one proton and diquark of another pro- 
ton. 



IV. MULTIPLICITY DISTRIBUTIONS 
FOR DIFFERENT PROCESSES 

We suppose that hadrons multiplicity dis- 
tribution in gluon string is normal distribution 
(Fig. 3). 

Gluon pairs are sequentially produced in fi- 
nal state marked by dashed block. This produc- 
tion takes place until energy of pairs approaches 
hadrons mass. Then gluon pairs decays at several 
observable hadrons. 

Every diagram of this type gives certain 
number of secondary hadrons. Number of such 
diagrams is considerably large. There is no sup- 
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FIG. 4: Vacuum contributions to non single diffraction 
cross sections of pp and pp collisions 



pression on energy because exchange is of vector 
type. Contributions from this diagrams have the 
same order of magnitude. So secondary hadrons 
multiplicity as random variable has to obey nor- 
mal distribution because of central limit theorem 
of probability theory. Thus multiplicity distribu- 
tion in gluon string is normal distribution. 

We suppose that hadrons multiplicity distri- 
bution in quark string is negative binomial distri- 
bution (NBD). This is in good agreement with 
data on e + e~ annihilation at low energies where 
one quark string is produced with assurance. So 
multiplicity distribution in two quark strings is 
convolution of two NBD, namely NBD with dou- 
ble parameters, and in three quark strings it is 
convolution of three NBD, namely NBD with 
triple parameters. 

Experimental charged multiplicities are nor- 
malized on non single diffraction cross sections, 
<?nsd = &tot — &el — &sd- Pomeron contributions 
are the same as for total cross sections 



nsd 



a,, 



CT o 



nsd ( 1 i rn 
(1 + d 1 



nsd i i rnsd i 2 \ 

ms + o 2 m s), 



here we divided formulae by constant term for 
convenience. 

The fitting of vacuum contributions to non 
single diffraction cross sections is shown on Fig. 4, 
data were taken from 0, E|. The value of x 2 t° 
number of degrees of freedom 14 to 8 is caused 
by difference between data of UA4 and UA5 col- 
laborations. 



0.011 and 5% sd 



The obtained coefficients S nsd = 0.075 ± 
0.007 ± 0.001 determine the 
weights of configurations with gluon string, two 
quark strings, three quark strings. 

Weight of normal distribution is defined by 
formulae and it is the same for pp and pp 
collisions. 



(3) 



l + $nsd \ ns + s nsd (\ns) 2 



Weights of double NBD are different for 
pp © and pp © interactions because in case of 
proton-antiproton configuration with two gluons 
also gives contribution to three quark strings. It 
is taken into account by coefficient c, which does 
not depend on energy of collision. 



8^ sd In s + 6^ sd (Ins) 2 
I + 5^ sd In s + 5^ sd {In 



8 nsd lns + {l-c)5^ sd (lns 



(4) 



(5) 



Therefore weight of triple NBD in configu- 
ration with three quark strings in pp collision is 
defined by formulae ©. 



c5% sd {\ns) 



1 + 5^ sd lns + 5% sd {lns) 2 



(6) 



But coefficient c can not be evaluated from the 
existing data, so we consider two possible values 
of it - 0.25 and 0.75. 



V. EXPERIMENTAL DATA ON 
CHARGED MULTIPLICITY 
DISTRIBUTIONS 

With assumptions stated in previous sec- 
tions we described experimental data on charged 
particle multiplicity for pp collisions at energy 
range from yfs = 44.5 to 200 GeV [5L la] and for 
pp collisions from 200 to 1800 GeV - 0. 

In case of pp interaction weights of distribu- 
tions are strictly defined by formulaes (|3|) , (j3]) . In 
Fig. 5 multiplicity distributions in pp interaction 
at yfs = 44.5 and = 62.2 GeV are shown. 
Here dash-dot line presents normal distribution 
in gluon string and dotted line presents double 
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FIG. 5: pp interaction at y/s = 44.5 {x 2 /ndf = 13/15) 
and vs = 62.2 GeV (x 2 /ndf = 13/16) 
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FIG. 6: pp interaction at y/s = 200 GeV 

NBD in two quark strings, solid line is sum of 
these distributions. 

But in case of pp interaction we have ad- 
ditional configuration with three quark strings, 
so weights of distributions are defined by formu- 
laes ©, (JSJ and The last two formulaes 

have undefined parameter - coefficient c which 
describes appearance of three quark strings. 

In Fig. 6 two fittings of pp collision at y/s = 
200 GeV are shown. Here dash-dot line presents 
normal distribution in gluon string, dotted line 
presents double NBD in two quark strings and 
dashed line presents triple NBD in three quark 
strings, solid line is sum of these distributions. 



FIG. 7: pp interaction at yfs = 1800 GeV 
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FIG. 8: pp interaction at y/s = 900 GeV (points) and 
prediction for pp at the same energy in LCNM for 
c = 0.75 (line) in logarithmic and linear scale 

On the left side the case when 25 per cent of two 
gluons give three quark strings is shown, on the 
right side - the case when 75 per cent of two glu- 
ons give three quark strings. Values of x 2 to num- 
ber of degrees of freedom are the same for both 
cases, x 2 / n df = 4/29. As energy increases higher 
value of coefficient c gives slightly better value of 
X 2 /ndf, for example, at the highest available en- 
ergy yjl = 1800 GeV x 2 /ndf = 138/115 when 
c = 0.25 and x 2 /ndf = 126/115 when c = 0.75 
(Fig. 7). 

After fitting all data we obtained energy 
dependence of parameters of gluon and quark 
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FIG. 9: Prediction of multiplicity distribution for pp 
interaction at y/s = 14 TeV in LCNM for c = 0.75 
(solid line) and for c = 0.25 (dashed line) in logarith- 
mic and linear scale 

strings ad thus we can predict the multiplicity 
distribution for pp collisions at 900 GeV. This is 
the first high enough energy for which it will be 
possible to compare multiplicity distributions for 
pp and pp interactions directly after start up of 
LHC. In the case of c = 0.25 the shapes for both 
reactions are practically the same. But in the case 
of c = 0.75 the systematic difference may be ob- 
served in region of "shoulder" given very precise 
measurements (Fig. 8). 

VI. PREDICTIONS FOR LHC ENERGY 
14 TEV 

We predict the value of total cross section at 
y/s = 14 TeV to be 

a vv ot = 101.30 ± 6.65. 

The multiplicity distribution for pp interac- 



tion can be described as sum of normal distri- 
bution in gluon string and double negative bino- 
mial distribution in two quark strings. Weights 
of these distributions are strictly defined by for- 
mulaes ([3j), Parameters of normal distri- 

bution are mean multiplicity in gluon string n g \ 
and standard deviation o g \ . Parameters of double 
NBD are mean multiplicity in quark string n q and 
shape parameter k q . Since parameters of quark 
and gluon strings must be obtained at simultane- 
ous fitting of pp and pp distributions, they con- 
tain uncertainty of three quark strings coefficient 
c. So we have two sets of energy dependencies for 
gluon and quark string parameters. 

c = 0.25 

n q = 3.45 + 0.46 In 2 y/s 
fc-i = -0.24 + 0.11 In y/s 
hgi = -5.86 + 3.77 In yfs 
<jgi = -2.67+ 1.52 In V* 
c = 0.75 

n q = 3.73 + 0.41 In 2 y/s 
k^ 1 = -0.30 + 0.13 In y/~s 
n gl = -6.65 + 3.97 In y/7s 
(jgi = -2.75 + 1.55 In y/s 

Predictions for charged multiplicity distri- 
bution for energy y/s = 14 TeV are shown in 
Fig. 9 for both cases of three quark strings co- 
efficient c. The values of mean multiplicity are 
the following 

/ 68.59 + 4.47, if c = 0.25; 
n PP ~ \ 63.20 + 4.19, if c = 0.75. 

It will be possible to distinguish between two 
values of coefficient c when experimental data are 
available. 
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